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We present a simple quantum theory for the pulsed light generated by a synchronously pumped 
optical parametric oscillator (SPOPO) in the degenerate case where the signal and idler trains of 
pulses coincide, below threshold and neglecting all dispersion effects. Our main goal is to precise 
in the obtained quantum effects, which ones are identical to the c.w. case and which ones are 
specific to the SPOPO. We demonstrate in particular that the temporal correlations have interesting 
peculiarities: the quantum fluctuations at different times within the same pulse turn out to be totally 
not correlated, whereas they are correlated between nearby pulses at times that are placed in the 
same position relative to the center of the pulses. The number of significantly correlated pulses 
is of the order of cavity finesse. We show also that there is perfect squeezing at noise frequencies 
multiple of the pulse repetition frequency when one approaches the threshold from below on the 
signal field quadrature measured by a balanced homodyne detection with a local oscillator of very 
short duration compared to the SPOPO pulse length. 
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T^ 1 , I. INTRODUCTION 

a: 

-t— > ■ Optical parametric oscillators (OPOs) are well-known and efficient sources of non-classical light. Whereas c.w. 

OPOs have been extensively studied both by theoreticians and experimentalists [l[, the quantum properties of pulsed 

OPOs have been so far much less investigated, in spite of their interest [2(- In particular, synchronously pumped 

rr-i OPOs (SPOPOs), iii which the pump pulses are temporally separated by the round trip time of the OPO cavity, 

>— — i, seem particularly promising, as theoretically shown in[3|, |j]. In these devices, the efficiency in twin-photon generation 
is enhanced twice, because of the high peak power in light pulses and because of field enhancement in a resonant 
cavity. In addition, the light emitted by SPOPOs has been theoretically shown to be either multi-mode squeezed 

_i_ ■ or multipartite entangled, which makes it an interesting resource for parallel transfer and processing of quantum 

,_^ | information |5( and for quantum metrology in the time domain, for example to measure ultra-short time delays with 

l/~j . a very high sensitivity [6[ . 

A detailed quantum analysis of the SPOPO has been performed in Refs. [j, |7|. In these papers the intracavity field 

r — . is considered in the frequency domain, i.e. expanded as a combination of longitudinal modes of the SPOPO cavity. It 
was demonstrated that a SPOPO emits a tensor product of squeezed "super-modes", each being a particular coherent 
superpositions of longitudinal modes of different frequencies. In this article we will take a rather different approach 

,- ^ ■ of the same problem, based on a description of the device in the time domain instead of the frequency domain. It is 
well-known that a theoretical description of pulsed processes is often more pictorial in such a frame, and the results 
obtained are complementary from the ones deduced from the frequency approach. More precisely, we will use the 
two-time technique used in Ref. |8| for the pulse laser generation. 

The paper is organized as follows. In Sec. [II] we present a physical model of the SPOPO in terms of two coupled 
Heisenbcrg-Langevin equations that describe in the time-domain the SPOPO operation both below and above oscilla- 
tion threshold. In Scc. lIIII wc determine the SPOPO oscillation threshold. In Sec. II VI we construct the pair correlation 
functions for the signal quadrature components outside the cavity in the below threshold configuration. In Sec. |V] 
the current correlation functions and their spectra are considered with and without time averaging by the detection 
process. Some interesting quantum features of the parametric oscillation are considered in the spectral domain in 
Sec. [VI] and discussed in Sec. IVIII In the Appendices a way to construct the main equations is considered in detail 
and typical physical parameters are selected. 



II. INTRACAVITY TWO-TIME DESCRIPTION OF FIELD PULSES 

Fig. 1 presents a sketch of a SPOPO that we consider here. A x^ nonlinear parametric crystal is inserted in an 
optical ring cavity and is pumped by a train of laser pulses with mean frequency u> p . We assume that the duration of 
the pump pulses t p is much shorter than the cavity round trip time Tr (under real conditions t p is about 10 -5 Tr). 
In the crystal takes place a degenerate type I parametric conversion of the pump field into a collinearly propagating 
signal field with mean frequency w p /2, as well as the reverse process. We assume a degenerate parametric interaction 
for the carrier frequencies of the pump and signal modes, meaning that the following phase matching condition is 
fulfilled: 

Ak = k p (uj p )-2k s (uj p /2) = (1) 

where k p and k s are wave vectors of carriers of pump and signal fields. 

The crystal is supposed to be thin so that we can neglect the influence of mismatch and dispersion of the group 
velocities on parametric interaction. Changes of fields caused by the parametric interaction are considered as pertur- 
bations of free propagation inside the crystal. In Appendix [X] we consider the physical parameters of the different 
parametric crystals and demonstrate that, for thin enough crystals, realistic conditions exist where we can avoid the 
undesirable influence of the dispersion phenomenon. 

We suppose that the SPOPO ring cavity is resonant and of high-finesse both for pump and signal fields, so that the 
doubly-resonant configuration is realized. We assume also that the cavity is dispersion-compensated by intracavity 
dispersive elements. This implies that optical pulses of arbitrary shapes are not distorted after one round trip inside 
the cavity, and also that it takes the same time Tr for pump and signal pulses to make a single round trip inside 
the cavity. In addition, the device is synchronously pumped, meaning that this time is equal to the pump repetition 
rate: Tr = T p . These hypotheses are set for simplicity, and the equations that we use in this paper could be easily 
modified in order to take into account deviations from these ideal conditions. 

In the degenerate parametric generation configuration that we consider here the field operator inside the cavity is 
equal to: 

E(z,t) = E p (z,t) + E s (z,t). (2) 

We use the plane wave approximation, so that the field amplitudes depend only on one longitudinal coordinate z 
measured along optical axis of the cavity. As the pump of the parametric crystal is realized by a train of optical 
pulses of duration close to 100/s it is possible to disjoint in the standard form quick oscillations of fields with optical 
frequencies uj p . s from slow changes of their envelopes [9[. At the crystal entrance for the pump p and signal s waves 
the two fields read 

ft \ 1 I 2 
Er{z,t) = i(-^—) e-^ t+ik ^A r (z,t), r = a,p, (co p = 2oj s ), (3) 

where n r = n r {ui r ) are the indices of refraction of the crystal and k r = k r {uj r ) are the pump and signal wave vectors 
determined by linear dispersion of the crystal. The slow amplitudes A r are normalized so that the mean values 

(A r A r ) have the meaning of mean fluxes in photons per second for the light beams of area S. One supposes that 
the crystal is placed just after the coupling mirror and that the field at the input boundary of the crystal (at z = 0), 
taking into account the periodic structure of the fields, is 

A r (0,t)=J2^r,n(-t-nT R ). (4) 

n 

Here A r _ n (t — tiTr) is the envelope of the n-th pulse. When t — jiTr appears as argument of the envelope we will treat 
it as time deviation from center of the pulse, i.e., it changes in the interval [— Tr/2, +Tr/2]. 

In order to describe SPOPO operation inside the cavity we use the two-time approach applied by Haus in Ref. Q 
for developing a quantum theory of actively mode-locked lasers. We assume that the envelopes of pulses are not 
significantly changed from one pulse to the next, a hypothesis that is typically valid in experiments with high-finesse 
cavity and weak parametric amplification. Then the dependence on discrete number n could be replaced approximately 
by a continuously varying temporal parameter T in the following way 

A rtn (t-nT R )^A r {t,T). (5) 



In the Appendix [B] we show in detail how to derive Hciscnbcrg-Langevin equations for these envelopes. We show that 
in the thin crystal approximation and when one neglects the dispersion phenomenon these equations read 



dA p (t,T) 
dT 



= -k p (A p (t,T) - A (t)) - gA s 2 (t,T) + F p (t,T), (6) 

^f^ = ~ K Mt,T) + 2gA p (t,T)Al(t,T) + F s (t,T) (7) 

Here k p and k s are the loss rates of pump and signal fields respectively; g is a constant characterizing the parametric 
coupling; Ao(t) is the classical steady-state envelope of the pump pulses inside the cavity, depending only on time t 
since the pump pulses are supposed to be perfectly identical. 

The Langevin noise sources F p (t,T) and F s (t,T) are characterized by the nonzero pair correlation functions: 

(F r (t,T)F}(t',T'))=2K r 5(T-T')6(t-t'), r = s,p. (8) 

III. OSCILLATION THRESHOLD OF THE SPOPO 

In order to describe the parametric effect below threshold, we can neglect pump depletion, and consider only the 
equation for the signal field, while the pump field remains in the coherent state with c- number amplitude Ao(t) 
imposed by the external pumping of the SPOPO. According to our model there is no phase modulation of the pump 
field, therefore without loss of generality we will treat the envelope of the pump pulse as a real positive quantity 
Aq = Aq > 0. Let us rewrite Eq. (J7J using the quadrature components of signal field 

X 8 (t,T) = \ (A s (t,T) + At(t,T)) , Y s (t,T) = I (A s (t,T)-Ai(t,T)) . (9) 

One then obtains two independent equations for the X s - and the Kj-quadraturcs 



d_ ( X s (t,T) \ / X s (t,T) \ / F>(t,T) 

dT { Y s (t,T) J- T[> { Y 3 (t,T) ) + { F>'(t,T) 



Mt) • + Sff^), do) 



where 



F s (t,T)=F&,T) + iF?(t,T), K± = K a (l±ii(t)), H(t) = — Mt)- (11) 

The coefficients k_ and k + are here time dependent, and can be negative in some time interval, leading to a divergence 
of the field quadratures, i.e. to a bifurcation of the system to the generation of "bright" pulses. The requirement 
k_ = determines therefore the threshold of the parametric oscillation. In order to stay below threshold we require: 

N (t)=A*(t)<N th = ^. (12) 

It is important to note that for the SPOPO the threshold depends not on the mean power of the pump as for the 
continuous generation but on the peak power of the pulse, because of the instantaneous character of parametric 
interaction. This is interesting for experimentalists, because the corresponding mean power at threshold can be 
extremely small: 

i r+Tn/2 
N = — dtN (t) « N th (13) 

J-R J-T R /2 

IV. BELOW THRESHOLD OPERATION: CORRELATION FUNCTIONS FOR THE FIELD 

QUADRATURES OUTSIDE THE CAVITY 

Let us now determine the quantum fluctuations and correlations of the output signal field quadratures, defined by: 

A° s ut {t) = X° ut (t) + iY° ut {t) (14) 



They are coupled to the intracavity amplitude A s (t) by the boundary condition on the output mirror, which has the 
form 

AT\t) = VT S A s (t) - ^K s A vac (t), (15) 

where T s , 1Z S (J~ s + 7Z- S = 1) are the reflection and transmission coefficients of the output mirror of the cavity (with 

n s « i). 



The formal solutions of Eqs. (1101) read 

W) \ = f T / F>(t,T>) \ e -^(T-T'). (16) 

Y s (t,T) J 1 M \F?(t,T')J { ' 

The two-time representation was introduced only as an intermediate operation to obtain the expression of the intra- 
cavity pulses by solving a differential equation. We can now come back to the description of the successive pulses 
using: 

T^nT R , t^t-nT R , X s (t,T) -> X s , n (t - nT R ), Y 8 (t,T) ->£,„(*- nT R ). (17) 

We thereby get the intracavity quadratures which have to be used in Eq. (|15[) 

X s , n (t - nT R ) \ f lTR , ( F' 8 {t,T>) \ - KT (nT R - T') 



dT > I - S V>- I \ -r, T yni R - ± ) , 18] 

Y s , n (t-nT R )J J^ \F> s <(t,T)) e 

These expressions enable us to determine the correlation functions of the output field quadratures in terms of simple 
temporal integrals: 

(X°%(t-nT R ) X°$(t' -n'T R )) 

(Y°£(t-nT R )Y°%(t'-n'T R )) 

= - A 5 nn ,5{t - t') ± Ks T R ptl^ML e -^(t) T R \n - n'\ 5{t _ t >_ {n _ „/ )Tfl) . (19) 

4 l^/j,{t- nT R ) 

The first term is due to incoming vacuum field reflected from coupling mirror of the cavity. The second term is related 
to the signal field coming out of the cavity. We see that at discrete times t ~ nT R the X-quadrature variance of a 
given pulse (when n = n') is increased above the vacuum level, and the Y-quadrature is squeezed. The presence of the 
delta-function in the expression is a sign that the different temporal parts of an individual pulse are not correlated, 
i.e., that the stretching/squeezing observed at different times are independent. The importance of the effect is defined 
by instantaneous value of the pumping amplitude, i.e., by the pump parameter /x(i). It is actually small because it is 
proportional to magnitude k s T r <C 1. 

Let us note that our model predicts quantum correlations between different pulses (when n =/= n') for the X- 
quadrature of the field, and anticorrelations for the Y-quadrature. The number of significantly correlated successive 
pulses can be evaluated by the factor in the exponential, which is proportional to {n s T R )~ l and roughly equal to 
cavity finesse at signal frequency. This result has the following simple interpretation [10(: the pump photons are 
paramctrically down-converted into pairs of correlated signal photons. The photons of each pair may leave the cavity 
in different pulses, during an overall time of the order of kJ 1 , giving rise to temporal correlations on the same range 
of time difference. The delta- function shows that the correlations between different pulses have a "local" character: 
they are effective only when the time differences are a multiple of the period T R . The equal time correlation of the 
signal field is obviously a consequence of the thin nonlinear crystal approximation used in the work. 

To measure such quantum effects on the field quadratures, one must use a balanced homodyne detection technique, 
which will be considered in the next section. 

V. BALANCED HOMODYNE DETECTION IN THE TIME DOMAIN 

Let us now investigate the measurement of field quadratures, obtained by a balanced homodyne detection of the 
output signal field (see Fig. 2). In this case, the current operator has the well-known form 

i(t) = 2(3(t) [cos$ X° u \t) + sin* F s outt (t)l , (20) 



where /3(t)e s * is the complex amplitude of the local oscillator. Choosing $ = and $ = tt/2, one can follow both 
quadrature amplitudes in the form 



l(t)=2j2/3(tn) 



Y out (t ) 

1 s,n \°nj 



t„ =t-nT, 



R 



(21) 



We supposed here that the local oscillator pulses arc all identical, that they do not have phase modulation and that 
their period is equal to the period of signal pulses T R . This implies that the envelope does not depend on pulse 
number n. We have made the same assumptions for the SPOPO pumping field. The parameters that can be changed 
in the homodyne detection setup are the duration of the local oscillator pulses tlo as well as their delay At relative 
to the signal pulses. 

Substituting here the expression (|19p . we derive the pair correlation functions for the currents in the most general 
form 

cm i(?))x 

fat) l{t'))y 

Y^f3 2 (t-nT R ) 



S(t - 1>) ± Ks T R 1 2 "(* - nT f . e-ttrfa) I* - *'l V 6 (t - t' - nT R + n'T R ) 
1 =F U(t - nTn) ^ 



(22) 



Here as before we can conclude that in the time domain the quantum effects are of the order k s T r and are therefore 
very small. We will see in the next section that they appear more clearly on the noise spectra. 

A real detector has a finite response time Td, and time averages the photodetection signal, so that the observed 
photocurrent I{t) is given by 



t+T D /2 

f(*) = -r / dti(t). 

t-T D /2 



(23) 



Eq. (|22p is valid when Td is much less than T R . Let us now consider the case, important in practice, when 
T R ^C Td ^C k^T 1 . This means that the pulse structure peculiar to the SPOPO is averaged. We expect therefore 
results looking like the c.w. regime. One gets in this case: 



where 



+Tr/2 



{I) = T- R 



(Ht) Ht'))x 
(Ht) Hnh 

dTp 2 (r), A T (t) 



-T R /2 



(I)S(t-t')±K s A T (t-t'), 



+Tr/2 
1 [ A (fit ^J^lL ' 

t~ r J dTp ^T^K7f 

-Tr/2 



(24) 



■ Ks (lTu(T))\t\ 



(25) 



By comparison with Eq. (|22p . wee see that only the equal time feature survives. Let us now choose the pulse of the 
local oscillator such that ft as a function of r is much narrower than /tx(r); then it is possible to use the approximation 
/3 2 (t) = Po <Kt), so that finally: 



(Ht) i(t'))x 



= (i) 



5(t-t')± Ks 2 ^°) e ^(lTM(0))|t-n 



(I(t) I{t'))y 

Close to threshold (1 - /i(0) < 1) we have: 

(i(t)i(t')) x /{i) = s(t-t') + 



1 T A*(0) 



(I) = f%. 

K ' T R 



1- M (0) 



,(l-/ i (0))|t-f| 



(I(t)i(t')) Y /(I)=S(t-t')-K s e 



-2K.\t-t'\ 



(26) 

(27) 
(28) 



VI. QUANTUM FEATURES ON NOISE SPECTRA 

In this section we determine the frequency spectrum of the quantum noise on the y-quadrature. This spectrum is 
defined as: 

+T/2 +T/2 

(i 2 ) fi = T lim i J dt J df <J(t)i(f)>ve in (*-*0. (29) 

-T/2 -T/2 

In the case of To ^C Tr, and after substituting Eq. (1221) into ((29)) . the spectrum reads: 

(z 2 ) =— / di/3 2 (i)(l- V P^ J. (30) 

V '° r W V m= fe 2i ...^(l + M*)) 2 + (^-27rm/T R )V 

This expression is valid for arbitrary shapes of pump and local oscillator pulses [i(t) and /3(t). 

In the case of very short local oscillator pulses probing signal ones at their peaks we have the simplified expression: 

(* 2 ww = i- E 42 '; /TM 2 - (3D 

m=o,i,2,- 4k s + (° - 27rm/r H ) 

One can see that the quantum noise reduction takes place not only in the vicinity of zero frequency, but also around 
all resonant frequencies of the cavity 2irm/Tn, as known also for the c.w. regime: indeed in the first experiment 
on generation of squeezed light based on four- wave mixing inside optical cavities |2| the measured photon pairs were 
symmetrically shifted with respect to the pump frequency by three cavity mode-spacing frequencies. In [111 ] and 
[12| squeezing at multiple longitudinal modes of a c.w. degenerate OPO was investigated theoretically and observed 
cxp er iment ally. 

The actual frequency response of detectors can be taken into account by multiplying the spectrum (|3Tj) by the 
detector frequency response function [llf . Obviously a high bandwidth photodetection set-up must be used to observe 
the noise reduction around multiple resonant frequencies. However if Tr <C Te> < kJ 1 , then in the spectrum only 
one quantum feature survives around the zero frequency and its spectrum is given by the well-known formula 

( ,! w"^-iip- (32) 

Turning back to the general expression (|30[) . one sees that it is nothing else than the time averaged noise spectrum 
of the signal field quadrature of a c.w. OPO pumped by a c.w. field of instantaneous value fJ.(t), weighted by the 
intensity of the local oscillator pulse. This is a consequence of the fact that different temporal parts of individual pulses 
of SPOPO are not correlated in our simple model and their nonclassical properties are defined by the instantaneous 
pump power. Consequently, the detected quantum noise is sensitive to the temporal properties of the local oscillator 
pulses, particularly to their duration and to the delay relative to signal pulses (see Figures [3] (a) and (b) respectively). 
As expected the model predicts maximal noise reduction for very short local oscillator pulses ideally synchronized 
with the signal maxima. 

VII. CONCLUSION 

In this paper, we have investigated the non-classical properties of a SPOPO operated below threshold using a time- 
domain approach, in the simplified case when it is possible to neglect the dispersion phenomenon. We have established 
differential equations for the intracavity field operator amplitude and on this basis investigated the quantum statistical 
properties of the signal field. As expected the typical region of correlation is the same as for the c.w. regime. However 
there arc supplementary details related to the pulse field structure, namely the fact that different times inside the 
limits of the same pulse turn out to be uncorrelated, whereas small but nonzero inter-pulse quantum correlations 
occur between fluctuations at times similarly placed in the different pulses. 

We have shown that the " instantaneous" homodyne signal (i.e. using a local oscillator of very short duration) turns 
out to be perfectly squeezed when approaching the oscillation threshold from below and at noise frequencies 27rm/Tjj 
(m = 0, 1,2, • • • ), just like in the c.w. case. From a practical point of view the major advantage of the SPOPO is its 
very low oscillation threshold in terms of mean pump power. This implies that one can use a moderately resonant 



cavity, with a high escape efficiency, and still have a threshold that can be reached using available mode-locked lasers 
with average powers in the lOOmW range. 

Let us mention finally that the Heisenbcrg-Langevin equations that we have established can be modified straight- 
forwardly to take into account experimental effects such as phase modulation and carrier-envelope phase shift of pump 
pulses, cavity detuning, and singly resonant operation. 
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Appendix A: Some quantitative estimations 

In this appendix we estimate the effects that occur in the propagation of femtosecond pump and signal pulses 
through a dispersive \ nonlinear crystal. These processes are the reshaping and chirping of pulses due to group 
velocity dispersion of the crystal and the pulse temporal walk-off due to mismatch of group velocities. The efficiency 
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of each process can be described by a characteristic propagation distance [14( . The group velocity mismatch of the 
pump and signal fields (v p and v s respectively) is characterized by the propagation distance L GV necessary for a 
signal pulse to walk through the pump one 

L GV = T P (A1) 

\i/v a -i/v p \' k ; 

where t p is the duration of pump pulses. The linear dispersion of the crystal is characterized by the distance L D over 
which the initial duration of Gaussian pulse t P:S will increase by a factor of V2 

L °--k. m 



NL 



The efficiency of parametric interaction is characterized by the distance L 

L NL = ^-r- (A3) 

2<r\A p eak \ 

where A^ ea - pump peak amplitude defined by external pumping as well as duration of pump pulses. For straight- 
forward parametric amplification with a pump peak amplitude the signal wave amplitude grow by a factor of e over 
this propagation distance. 



In Table 1 we present characteristic distances evaluated for the following experimental parameters considered in [7| : 
pump wavelength A p = 0.4/^m, duration of pump pulses - r p = 100/s, c.w. threshold pump power inside the cavity 
- P t h = 25mW; crystal length - I = 0.1mm. We have also supposed following duration of signal pulse r s = 100/s in 
order to estimate Lf using (JA2I) . 



crystal 


t GV 

L , mm 


L®, mm 


L®, mm (r s = 100 fs) 


L NL , mm 


BBO 


~ (0.25 -Hi) 


- (10 -=-20) 


-(20^30) 


-4.7 


KNb0 3 


~ (0.1 -=-0.75) 


- (10 -=-50) 


- (20 -=- 40) 


- 1.6 



Table 1: Estimations of characteristic distances 
Thus data of the Table 1 show that for used experimental parameters the following relation holds 

I < L GV < L NL < L D 



(A4) 



Therefore, in the cases we consider here all the discussed processes can be considered as small perturbations of free 
propagation of pump and signal fields inside the crystal (both BBO and KNbOs). More precisely the dispersive 
reshaping and chirping of envelopes of fields are the weakest processes. 

Expression (|A2[) also defines characteristic time for the crystal of given length I 



vw 



(A5) 



Thus, neglecting crystal dispersion we can correctly describe the field properties only at time scales larger than t® . 
For the present experimental parameters this time is on the order of 6 fs. 

Appendix B: Derivation of the time-domain equations for the intracavity pulsed field 

Let us write the pump and signal fields propagating through the crystal in terms of slowly varying complex envelopes 
A r (z,t) introduced by means of the definition 



E r (z, i) = i 



hu r 



2n r eocS 



1/2 



—i(jj r t-\-ik r z 



A r (z,t), r = s,p, (lo p = 2lo s ). 



(Bl) 



Here k r = k r {u> r ) are the pump and signal wave vectors determined by linear dispersion of the crystal; n r = n r {uj r ) - 
the indexes of refraction of the crystal. 

Collincar propagation of pump and signal fields inside x nonlinear crystal with dispersion is described by following 
two coupled traveling-wave equations written for operators of complex envelopes of fields [14| 



-iAkz 






Here a is the coupling constant, proportional to the nonlinear susceptibility of the crystal x 



a 



r (2) / tkv p ui 2 s 



c y 2eQcSn p n 2 s 
The inverse group velocities and the dispersions of group velocities are respectively 

-l2, 



_! _ dk r (uj) 



Pr 



d'k r (uj) 



(B2) 

(B3) 

(B4) 
(B5) 



Exponents on the right hand sides of equations describe phase mismatch of pump and signal carrying waves along 
their propagation inside the crystal 



Afc = fcp(wp) - 2k s (u) s ). 



(B6) 



In this paper we consider degenerate regime of SPOPO operation, Ak = 0. 

These equations assume an instantaneous nonlinear response of the medium and use the approximation of slowly 
varying envelopes. Also linear dispersion of the crystal is considered up to the second order in the perturbative 
expansion around the pump and signal carrying frequencies |15| . Under these approximations equations do not 
describe arbitrarily fast dynamics of envelopes of fields, but they are quite reasonable for pulses with duration of tens 
of femtoseconds. 

Let us assume that j3 r = 0, meaning that the parametric crystal is thin enough so that the dispersive reshaping 
and chirping of envelopes of fields along their propagation inside the crystal are negligible. Under this approximation 
second-order time derivatives in equations (|B2[) and (|B3|) are dropped. We also assume the parametric interaction is 
very weak in a single pass through the nonlinear crystal, so that we can use in the right hand sides of the equations 
the amplitudes of the free propagating fields. In addition, we will assume for simplicity that the group velocities of 
the pump and signal fields are equal: v p = v s = v. The traveling-wave equations take now the simplified form: 

-^ + ~)A p (z,t) = -aA 2 s (0,t-z/v), (B7) 



TTz + Idi) ^'^ = 2aA p {0,t-z/v)At(0,t-z/v) 



(B8) 



The solutions of these equations can be easily obtained by replacing the time t by its reduced value rj = t — z/v and 
integrating obtained equations over the crystal length. The solutions read 

A p {l, t) = i p (0, t - l/v) - al A 2 S {0, t - l/v), (B9) 

A s (7, t) = A s (0, t - l/v) + 2al A p (0, t - l/v)A\ (0, t - l/v) (BIO) 

Now let this thin parametric crystal be placed inside the high-finesse cavity just near the input mirror (see Fig. 1). 
Putting the coordinate axis z parallel to the optical axis of the cavity one can couple by the following relation the 
intracavity field amplitudes A r (L,t) and A r (0,t) before and after the coupling mirror respectively with the amplitude 
Aj. ra (0, t) of the external field entering the cavity 



A r (0,t) = y/K r A r (L,t) + y/% 4J?(0,t), (Bll) 

where lZ r and T r (r = p, s) arc reflection and transmission coefficients of the mirror. The slow amplitudes of the input 
fields defined in the same manner as intracavity ones obey following commutation relations Q 

[A i r n (z,t),Al n Hz,t')]=S(t-t'), [Al n (z,t),Al n (z,t')]=0. (B12) 

Using solutions (|B9[) - (|B10[) obtained under thin crystal approximation and assumption of equal group velocities it 
is straightforward to write following the relations that couple the amplitudes of fields before and after one round-trip 
inside the SPOPO cavity 



A p (L,t) = (A p (0,t-T R )-alA 2 8 (0,t-T R )j e iu " T ^, (B13) 

A s (L,t)= (A s (0,t-T R ) +2alA p (0,t-T R )Al(0,t-T R )) e iu ° T ^° (B14) 

Here T R - round-trip time of pump and signal pulses inside the cavity; T p h, p and T p ] ltS - round-trip times of pump 
and signal carrier waves. We assume that the pump and signal carriers are resonant for the cavity so that conditions 
^pTph.p = 27rm and w s T p h tS = 27m hold. 

Representing continuous slowly varying envelope of the signal field as the following piecewise function: 

l a (0 > t)=534,,„(t-nT«) (B15) 

and combining (|B11[) with (|B14|) one obtains 

A s , n (t - nT R ) = ^/n s (i s ,„_i(i - (n - l)T R ) + 2al i p , n _i(t - (n - 1)^)1+^ (* - (n - 1)Tr)) + 

+VT s A™ n (t-nT R ), (B16) 

where the time t is treated as the time deviation from the pulse center. We see that there are actually two time 
arguments, because the index n informs us about exchanges from pulse to pulse. Strictly speaking this second 



10 

argument is discrete with distance between adjoining point equal to Tr. However assuming that the pulse envelope 
does not change essentially after one round trip inside the cavity one can neglect the time interval Tr in comparison 
with other typical. This enables us to consider the time variable T as continuous. Such a two-time approach has been 
used in Ref. [8[ for description of pulsed lasers. We can then make the following replacements in Eq. (|B16|) 

if. t^autm \ A ^ f ~ nTfl ) ~ A s, n -i(t -(n- 1)T R ) dA a (t,T) 

A r , n {t-nT R ) -> A r (t,T) (r =p,s), > — (B17) 

Jfl ell 

We can finally use the same procedure to get the equation for the pump pulse amplitude. We have then justified in 
this appendix the introduction of equations (jB]) and (Q of the present paper. 



11 




Figure 1: Synchronously pumped OPO 
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Figure 2: Balanced homodyne detection 
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Figure 3: Quantum noise of the squeezed quadrature normalized to the shot noise at zero frequency as a function of duration 
of local oscillator (LO) pulses (a) and their delay relative to signal ones (b). Both times normalized to pump pulses duration. 
LO pulses and pump ones are Gaussian. At the Fig. (a) LO pulses are ideally synchronized with signal At = 0; at the Fig. (b) 
two fixed durations of the pulses are considered. In both cases n a Tn = 0.1 and fj,(0) = 0.9. 



